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Abstract 

We consider a finite element approximation of a general semi-linear stochastic partial differ¬ 
ential eqnation (SPDE) driven by space-time multiplicative and additive noise. We examine 
the full weak convergence rate of the exponential Euler scheme when the linear operator is self 
adjoint and provide preliminaries results toward the full weak convergence rate for non-self- 
adjoint linear operator. Key part of the proof does not rely on Malliavin calculus. Depending 
of the regularity of the noise and the initial solution, we found that in some cases the rate 
of weak convergence is twice the rate of the strong convergence. Our convergence rate is in 
agreement with some numerical results in two dimensions. 
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Strong convergence 


1. Introduction 

The weak numerical approximation of an Ito stochastic partial differential equation defined in 
D C is analyzed. Boundary conditions on the domain D are typically Neumann, Dirichlet or 
Robin conditions. More precisely, we consider in the abstract setting the following stochastic 
partial differential equation 

dX = {AX + F{X))dt + B{X)dW, X(0) = Xo, t e [0,T], T > 0 (1) 

on the Hilbert space L^(D). Here the linear operator A which is not necessarily selfadjoint, 
is the generator of an analytic semigroup S{t) := > 0. The functions F and B are 
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nonlinear functions of X and the noise term W{t) is a Q-Wiener process defined on a filtered 
probability space (D, X, P, that is white in time. The hltration is assumed to fulhll 

the usual conditions (see e.g. Dehnition 2.1.11]). For technical reasons more interest will 

be on a deterministic initial value Xq G H. The noise can be represented as a series in the 
eigenfunctions of the covariance operator Q given by 

W{x,t) = (2) 


where i E are the eigenvalues and eigenfunctions of the covariance operator Q 

and /3i are independent and identically distributed standard Brownian motions. Under some 
technical assumptions it is well known (see 0,0 01 ) that the unique mild solution of ([I]) is 
given by 


X{t) = S{t)Xo+ f S{t-s)F{X{s))ds+ f S{t- s)B{X{s))dW{s) 


( 3 ) 


Equations of type ([T]) arise in physics, biology and engineering [27|, |28|, l9j and in few cases, 


exact solutions exist. The study of numerical solutions of SPDEs is therefore an active research 
area and there is an extensive literature on numerical methods for SPDE ([T]) 0,0,0,0,0, 
00, 0 |. Basicaly there are two types of convergence. The strong convergence or pathwise 
convergence studies the pathwise convergence of the numerical solution to true solution while 
the weak convergence aims to approximate the law of the solution at a hxed time. In many 
applications, weak error is more relevant as interest are usualy based on some functions of the 
solution i.e, E$(X), where $ : if —)■ R and E is the expectation. Strong convergence rates 
for numerical approximations of stochastic evolution equations of type ([1]) with smooth and 
regular nonlinearities are well understood in the scientihc literature (see [ 0 , 14, 0 , 0 0 , 
34j . l28l. |22| | and references therein). Weak convergence rates for numerical approximations of 
equation (jl]) are far away from being well understood. For a linear SPDE with additive noise. 


the solution can be written explicitly and the weak error have been estimated in [2a, [2^ with 


implicit Euler method for time discretization. The space discretization have been performed 
with hnite difference method l 


24 and hnite element method 24 . The weak error of the 


implicit Euler method is more complicated for nonlinear equation of type ([T]) as the Malliavin 
calculus is usually used to handle the irregular term and the term involving the nonlinear 
operators F and B (see 0,0, i). In almost all the literature for weak error estimation, the 
linear operator A is assumed to be self adjoint. Furthermore no numerical simulations were 
made to sustain the theoretical results to the best of our knowledge. In this paper we consider 
a stochastic exponential scheme (called stochastic exponential Euler scheme) as in and 
provide the weak error of the full discrete scheme (Theorem [621 Theorem 16.31 and Remark l6.3ll 
where the space discretization is performed using hnite element, following closely the works in 
00 on another exponential integrator scheme. Our weak convergence proof does not use 
Malliavin calculus. Furthermore we provide some preliminaries results (weak convergence of 
the semi discrete scheme in Theorem 14.II and Theorem 15. ip toward the weak convergence when 
the linear operator A is not necessarily self-adjoint, and provide some numerical examples to 
sustain the theoretical results. Recent work in 0 is used to obtain optimal convergence order 
for additive noise when the linear operator is self adjoint in Theorem 16.31 and Remark 16.31 
We also extend in Theorem 16. II the strong optimal convergence rate provided in 1^, Theorem 
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1.1] to non-self-adjoint operator A. Note that as the operator A is not necessarily sef-adjoint, 
our scheme here are based on exponential matrix computation. The deterministic part of this 
scheme have be proved to be efficient and robust comparing to standard schemes in many 
applications [^, 28, 30, 29] where the exponential matrix functions have been computed using 
the Krylov subspace technique [l^ and fast Leja points technique [^. For convenience of 
presentation, we take A to be a second order operator as this simplihes the convergence proof. 
Our results can be extended to high order semi linear parabolic SPDE. 


The paper is organized as follows. Section |2] provides abstract setting and the well posedness 
of ([I]). The stochastic exponential Euler scheme along with weak error representation are 
provided in Section |21 The temporal weak convergence rate of the stochastic exponential 
Euler scheme is provided in Section 0] for additive noise and in Section 0] for multiplicative 
noise. Note that in this section the linear operator A is not necessarily self-adjoint. Section [6] 
provides strong optimal convergence rate of the semi discrete solution for non-self-adjoint 
operator A along with full weak convergence rate of the stochastic exponential Euler scheme 
for self-adjoint operator A. Numerical results to sustain some theoretical results are provided 
in Section [71 


2. The abstract setting and mild solution 

Let us start by presenting briefly the notation for the main function spaces and norms that 
we use in the paper. Let FT be a separable Hilbert space with the norm || • || associated to 
the inner product (,}h- For a Banach space 1/ we denote by || • ||n_the norm of U, L{U,H) 
the set of bounded linear mapping from U to H and by L2(]I3, f/jj the Hilbert space of all 
equivalence classes of square integrable [/—valued random variables. For ease of notation 
L{U, U) = L{U). Furthermore we denote by £i(f/, H) the set of nuclear operators from U to iF, 
C 2 {U, H) := HS{U, H) the space of Hilbert Schmidt functions from U to H and Cj^iU, H) the 
space of not neccessarily bounded mappings from U to H that have continuous and bounded 
Frechet derivatives up to order F, k E N. For simplicity we also write jCi{U, U) = C\{U) and 
C2{U,U)=C2{U). 

For a given orthonormal basis (e,) of U, the trace of I G jCi{U) is dehned by 

:='^{lei,ei)u, (4) 

while the norm of I E C 2 {U) is dehned by 

IKII£2(C/) ^ 

Note that the trace in(|3]) and the Hilbert Schmidt norm in (jS]) are independent of the basis 

(e*)- 


^ID) is the sample space 
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Let Q : H ^ H he an operator, we consider throughout this work the Q-Wiener process. We 
denote the space of Hilbert-Schmidt operators from to Lf by L 2 := H) = 

HH) and the corresponding norm ||.||i;,o by 


IK 



¥Q^^‘^\\c2{h) 



leLl 


Let ip ■. [0, T] xff —)■ L 2 be a valued predictable stochastic process with P ||(y9||^ods < 00 
1. Then Ito’s isometry (see e.g. j^, Step 2 in Section 2.3.2]) gives 


E|| / ^fdwf= / 


«e|0,r]. 


Let us recall the following proprieties which will be used in our errors estimation. 
Proposition 2.1. Q/ Let /, /i, I 2 he three operators in Banach spaces, the following proprieties 


hold 


• If I E Ci{U) then 


mi)\<\\i\\c,iu). (6) 

• If h E L{H) and I 2 E Ci{H), then both I 1 I 2 and I 2 I 1 belong to Ci{H) with 

Trihk) = Trikh). (7) 

• If h E € 2 ( 11 , 11 ) and I 2 E € 2 ( 1 !, U), then I 1 I 2 E €i(H) with 

WhkWciiH) < \\h\\c2{u,H)\\k\\c2{H,u)- ( 8 ) 

• If I E € 2 ( 11 , 11 ), then its adjoint I* E € 2 ( 11 , U) with 

\\^*\\c 2{H,U) = \\mC2{U,H)- (9) 

• If I E L(U,H) and Ij E €j(U), j = 1,2, then llj E €j(U,H) with 

\\llj\\cj(U,H) < |K||L(C/,H)|Kil|£,(C/), j = 1,2. (10) 


For classical well posedness, some assumptions are required both for the existence and unique¬ 
ness of the solution of equation ([T]). 

Assumption 2.1. The operator A : T>(A) <Z H ^ H is a negative generator of an analytic 
semigroup S(f) = t > 0. 

In the Banach space P((—A)"/^), a e M, we use the notation ||(—A)"/^.|| =: ||.||q,. We recall 
some basic properties of the semigroup S(t) generated by A. 
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Proposition 2.2. [Smoothing properties of the semigroup ]l^]] 

Let a > 0, /9 > 0 and 0 < 7 < 1, then there exist C > 0 such that 




< Ct ^ for t > 0 

< CL^ for t > 0 . 


In addition, 


{-AfS{t) = 
If (3 > 
\DlS{t)vy < 


S{t){-Af on V{{-Af) 

7 then V{{-Af) C V{{-Ay), 

^^-z-(/3-a)/2 t > 0, v e D((-A)“/2) / = 0,1, 


where D\ := —. 
* dt^ 


We describe now in detail the standard assumptions usually used on the nonlinear terms F,B 
and the noise W. 


Assumption 2.2. [Assumption on the drift term F[ There exists a positive constant 
L > 0 such that F : H ^ H satisfies the following Lipschitz condition 


\\F{Z)-F{Y)\\<L\\Z-Y\\ V Z,YeH. 


As a consequence, there exists a constant C > 0 such that 

||F(Z)|| < ||F(0)|| + ||F(Z)-F(0)||<||F(0)||+L||Z||<C(1 + ||Z||) Z e H. 


Assumption 2.3. [Assumption on the diffusion term B[ 

There exists a positive constant L > 0 such that the mapping B : H ^ II' 2 {.H) satisfies the 
following condition 


||B(Z) - B(Y)\U,(„) <L\\Z- y|| VZ, YeH. 


As a consequence, there exists a constant C > 0 such that 

\\B(Z)\\c,m < ||F(0)||£,(„| + ||F(Z)-B(0)|U.m 

< ||B(o)||£,(„, + L||Z|| <C(l + ||Z||) ZeH. (11) 

Theorem 2.1. [Existence and uniqueness (J^])] 

Assume that the initial solution Xq is an Fq— measurable H—valued random variable and 
Assumption \2.1[ Assumption \2.A Assumption 1^.31 are satisfied. There exists a mild solution 
X to ^ unigue, up to equivalence among the processes, satisfying 

||X(s)||2ds < cx)^ = 1. 
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( 12 ) 








For any p>2 there exists a constant C = C{p,T) > 0 such that 

sup E||X(t)f <C{1 + EllXof). (13) 

te[o,T] 

For any p > 2 there exists a constant Ci = Ci{p,T) > 0 such that 

Esup ||X(t)|r <C'i(l + E||Xo|n. (14) 

te[o,T] 

The following theorem proves a regularity result of the mild solution X of ([T]). 

Theorem 2.2. [Regularity of the mild solution (\2^])] 

Assume that Assumption \2. 1[ Assvmption \2.A and Assvmption \2. A hold. Let X he the mild so¬ 
lution o/ (C]) given in If Xq G L 2 (]D),'r’((—A)^/^)), ft G [0,1) then for all t G [0,T], X(t) G 
L 2 {B),V{{-A)h/^)) with 


(E||X(i)|l?)‘'Tc (l+(E||X„|||)‘'"). 


For the weak error represention, we will need the following lemma. 

Lemma 2.1 (Ito’s formula). Let (D, X, P; be a filtered probability space. Let 0 

and T he H—valued predictable processes, Bochner integrable on [0,T] P-almost surely (see 
0 /; , and Yq be an Fo-measurable H—valued random variable. Let G : [0,T] x P —)■ M and 
assume that the Frechet derivatives Gt{t,x), Gx{t,x), and Gxx{t,x) are uniformly continuous 
as functions of{t,x) on bounded subsets o/[0;T] xH. Note that, for fixed t, Gx{t,x) G L(P, M) 
and we consider Gxxit,x) as an element of L{H). Let W be the Q-Wiener process. If Y 
satisfies 

y(f)=y(0)+ f (j){s)ds+ f T(s)dlF(s), (15) 

Jo Jo 

then P-almost surely for all t E [0, T], 

G{t,Y{t)) = G{0,Y{0))+ fGx{s,Y{s))^{s)dW{s) 

t 

+ f {G,(s,F(s)) + G4s,y(sMs) + i2v(G,,(s,F(s))'r(s)Qi(<I>(s)Ql)-)} 

(16) 

A proof of this lemma can be found in j^. 


3. Application to the second order semi—linear parabolic SPDE 

We assume that hi has a smooth boundary or is a convex polygon of d = 1, 2, 3. In 
the sequel, for convenience of presentation, we take A to be a second order operator as this 
simplihes the convergence proof. The result can be extended to high order semi linear parabolic 
SPDE. 


ds 


6 





More precisely we take H = L‘^{Q) and consider the general second order semi-linear parabolic 
stochastic partial differential equation given by 

dX{t, x) = (V ■ DVX(t, x) - q ■ VX(t, x) + fix, X{t, x))) dt + b{x, X{t, x))dW{t, x), (17) 

X ^ fl,t E [0,T] where / ihlxM—)-Misa globally Lipschitz continuous function and 
feihlxM—)-Misa continuously differentiable function with globally bounded derivatives. 


3.1. The abstract setting for second order semi-linear parabolic SPDE 

In the abstract form given in CQi. the nonlinear functions F : H ^ H and B -. H ^ 
HSiQ^/‘^iH),H) are dehned by 

iFiv))ix) = fix,vix)), iBiv)u)ix) = bix,vix)) ■ uix), (18) 

for all X G n, v E FI, u E Q^/^(i7), with H = L'^ifl). Note that we can also dehne 
B : H ^ C 2 iH) by 

iBiv)u)ix) = 6(x,n(x)) ■ Q^^^m(x), (19) 


for all X G n, v E H, u E H. 

In order to dehne rigorously the linear operator, let us set 


M = 


.A 

^ dxi I ^’^dxj) dxf 

i.j=l ^ ' 2=1 


( 20 ) 


where we assume that Dij E qi E and that there exists a positive constant 

Cl > 0 such that 


d 

'^Dijix)iiij > cil^l^ G X G FL, Cl > 0. 


( 21 ) 


We introduce two spaces El and V where El C Id depends on the choice of the boundary 
conditions for the SPDE. For Dirichlet boundary conditions we let 

1/ = e = B^(Q) = {vE Ff\Q) :v = 0 on 

and for Robin boundary conditions, Neumann boundary being a special case, we take Id = 
JF^(F) and 


See 


H = {x G 77^(17) : dvjdvji^ + aov = 0 on 972} , ao E 
for details. The corresponding bilinear form of —A is given by 

d 

'H~n 

OX 


aiu, v) = 


/A 

\io=i 


du dv 


dxj dxi 




u,v eV 


( 22 ) 
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for Dirichlet and Neumann boundary conditions, and by 


a{u, v) 



du dv 
dxj dxi 


d 

+ & 

i=l 



dx+ aou V dx, 


I an 


M, n e V, 


(23) 


for Robin boundary conditions. According to Garding’s inequality (see 
two positive constants Cq and Aq such that 


28 


there exist 


a(v,v) + cjlliif > Vti e V. 


(24) 


By adding and subtracting c^Xdt on the right hand side of ffT7|l . we have a new operator 
that we still call A corresponding to the new bilinear form that we still call a such that the 
following coercivity property holds 

a(v,n) > Ao||n||^i(n), Vn e R. (25) 

Note that the expression of the nonlinear term F has changed as we include the term —CqX 
in a new nonlinear term that we still denote by F. 

Note that a(,) is bounded in R x R, so the following operator A : R —)■ R* is well dehned by 
the Riez’s representation theorem 

a{u,v) =—{Au,v), \/u,v^V, (26) 

where V* is the adjoint space (or dual space) of R and (,) the duality pairing between V* 
and R. By identifying H to its adjoint space H*, we get the following continuous and dense 
inclusions 


V CH CV*. 


So, we have 


(27) 


{u,v)h = {u,v) Vm G if, Vf G R (28) 

The linear operator A in our abstract setting ([T]) is therefore dehned by fl26p . The domain of 
A denoted by T>(A) is dehned by 

V{A) = {m G R, Am G H}. (29) 

We write the restriction of A : R ^ R* to T>{A) again Iw A, which is therefore regarded as 
an operator of H (more precisely the if realization of A |8|, p. 812]). The coercivity property 
implies that A is a sectorial on L^(G) i.e. there exists Ci, 6 G (|vr, tt) such that 


ll(-^-^“"4) < -pi A G fe. 


(30) 


where fe = {A G C : A = p > 0, 0 < |(/)| < 
generator of bounded analytic semigroups S{t) : 


9'^ (see 11, [sl). Then A is the inhnitesimal 
= on F^iyi) such that 


Sit) := 


e*^(Af - A)-MA, 


1 

2m 


t > 0, 


(31) 





where C denotes a path that surrounds the spectrum of A. The condition property fl2S]) also 
implied that —A is a positive operator and its fractional powers is well dehned for any a > 0, 
by 


^a-ietA 


dt, 


(32) 


(-A)“ = ((-24)-i)-S 

where r(Q;) is the Gamma function (see 0 ). Functions in El satisfy the boundary conditions 
and with H in hand we can characterize the domain of the operator (—and have the 
following norm equivalence (0,00) for a G {1, 2} 


lTllH“(r2) 

p((_yl)«/ 2 ) 

V{-A) 


||(- 2 l)“/\|l =: ||n|U, Vn e 

H n (Dirichlet boundary conditions), 

H, V{{—Ay^‘^) = H^{Q) (Robin boundary conditions). 


3.2. Numerical scheme and weak error representation 


We consider the discretization of the spatial domain by a hnite element method. Let 7^ be a 
set of disjoint intervals of G (for d = 1 ), a triangulation of G (for d = 2 ) or a set of tetrahedra 


(for d = 3) with maximal length h satisfying the usual regularity assumptions la, llO 


Let 14 C R denote the space of continuous functions that are piecewise linear over the 
triangulation 71. Note that for high order polynomial, high order accuracy in space can be 
achieved. To discretize in space we introduce the projection Ph from L^(r2) onto 14 dehned 
for u G by 


{PhU, x) = {u,x), V X G 14- (33) 

The discrete operator Ah'.Vh^ 14 is dehned by 

(24h(p, x) = x) = x), X e 14. (34) 

Like the operator A, the discrete operator Ah is also the generator of an analytic semigroup 
Sh ■= Here we consider the following semi-discrete form of the problem ([T]), which 

consists to hud the process X^(t) = X^{.,t) G 14 such that for t G [0,T], 

dX^ = {AhX^ + PhF{X’^))dt + PhB{X^)PhdW, X'*(0) =: X^ = PhX^. (35) 

Note that fl35p is a hnite dimensional stochastic equation. The mild solution of (155]) at time 
tm = mAt is given by 

ptm 

X\tm) = Sh{tm)PhXo+ Sh{tm-s)PhF{X\s))ds 

Jo 

ptm 

+ Sh{tm-s)PhB{X\s))PhdW{s), (36) 

Jo 

■where At = T/M, m G {0,1,..., M}, M G N. 
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Then, given the mild solution at the time tm, we can construct the corresponding solution at 

tm+l tlS 


-At 


X^{tm+l) = Sh{At)X^{tm)+ Sh{At-s)PhF{X^{s + tm))ds 


Cm+1 


+ 


Sh{tm+i - s)PhB{X\s))PhdW{s). 


To build the numerical scheme, we use the following approximations [22 


P^,F{X\tm + s)) ^ PhF{X\tJ), se[0. At], 
Sh{tm+i-s)PhB{X\s)) ^ Sh{At)PhB{X\tJ), se[t™,Wi]. 

We can dehne our approximation X^^ of X{mAt) by 

Xt+i = ‘ - I) P,F(X!:J + )P» - W,,„) 

We introduce for brevity the notations 

Si(i)=e“*, «,;(() = Bk ■. X ^ PuB(X)Pu. 


(37) 


We can then rewrite the scheme fl?7|l as 

.Y*+, = SK(Xt)Xt + XtSl(At)P,F(Xt) + S4A()B„(At) - H')„) . (38) 

In order to study the weak convergence of the approximation of the solutions we dehne the 
functional 

/(t,^) = E [<l>(x"(t,^)], t e [0,T], ^ e W, (39) 

where X^{t,'ip) is dehned by fl36p with the initial value Xq = 'ip E 14- If can be shown (see 
Theorem 9.16 of [9]) that /i^(t,4) dehned by is diherentiable with respect to t and twice 
diherentiable with respect to ip, and is the unique strict solution of 

\ /r'^(0,4) = $(V'), ^ ’ 


We can here, using the Riesz representation theorem, identify the hrst order derivative of 
/i^(f, %p) with respect to ip, denoted as DjjAit, ip) with an element of 14 and the second deriva¬ 
tive denoted as D'^^^{t,ip) with a linear operator in 14- More precisely 

Dli^{t,ip){(pi) = {D^i^{t,ip),(pi)H, \/ip,(piEVh, (41) 

T>V^(t,4)(0i,02) = '0)01, 02)rr, V4, 01,02 e 14- (42) 


The following theorem, which is similar to 32|, Theorem 2.2] is fundamental for our convergence 
proofs. 
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Theorem 3.1 (Weak error representation formula for the semi discrete problem). 

Assume that all the conditions in the assumptions above are fulfilled and let {lT(t)}tg[o,r] be 
a cylindrical H—valued Q—Wiener process. Then for $ G M) the weak approximation 

error of the scheme in fl38l) has the representation 


E[>1>(A'S,)] - E[>l>(A"‘(r))] 

M -1 ( ptm+l 

= E / E 

m =0 KJtm 

/*^7n+l 


{Dfr\T - s, X\s)), PhFiXlf) - PhF{X\s))) 


H 


ds 


-iE 


Tr 


- s,.Y'-(s)) {S„(s - UB„{Xi,)y 

-Bft(A‘)B»(A*)-)] ds |, 


where X^{t) is a continuous extension of Xf^, defined by 

X\t) = Sh{t-tjX^ + {t-tJSl{t-tjPhF{Xi) 

+Sh{t-tjBf,{X^){W{t)-W{tJ), /orte 


(43) 

(44) 


Proof. Introduce the process : [0,T] x 14 — )■ M, given by 

fi) = p!^{t, 4fe(-t)4), (45) 

which is twice differentiable with respect to and by previous identifications satisfies 

L>z/^(t,4)0 = {Dp.’^{t,Sh{-t)fi),Sh{-t)(j))H, &Vh (46) 

D‘^iy^{t,fi){(pi,(j)2) = {D‘^ia^{t,Shi-t)fi){Shi-t)(j)i),Shi-t)(j)2)H, ^Vh. (47) 

One can then check that v^{t,f)) solves the equation 

^(t,4) = {Do\fi'ilj),Sh{t)PhF{Sh{-t)fi)H + \TT[D^p\t,fi)Sh{t)Bn{fi){Sh{t)Bhm% 
z/'*(o, 4 ) = <l)( 4 ), 4 e 14 . 

(48) 

Indeed, since = —AhSh{—t)fi, we have that 

^(t, 4) = ^(t, s„{-t)fi) + (Zi/i"(t, 44 - 44 ), - A44-t)4)H. 

Using the Kolmogorov equation fl4U]) and the fact that (, )h is symmetric yields 

^(44) = ( 4 , 44 - 44 + PhF(44-44),^/i"(4 44 -44))h 

+iTr [D^fi\t, 44-44)44-t)44-444t)44t)544)*] 

+{Dfi'^{t, 4,(-44), -AhSh{-t)fi)H- 

As D‘^fi^{t,'ijj) is identihed as linear operator in U,, using Proposition 12 . If mainlv relation (17|) ). 
the dehnition of the trace (jl]), the expression fl47|) and the fact that {,)h is symmetric allow 
to have 

^(t,4) = {Dfi^{t,Sh{-t)fi)Sh{-t)fi,Sh{t)PhF{Sh{-t)ilj))H 
+iTr [ziV4t,4)444S44)(444S44))*] 

= {Du'^it, 4), 44t)P,F(44-44)rf + iTr[PV4t, 4)44t)P44)(444544))*]. 
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Now let Z^(t) = ShiT — t)X^{t), using and the notation Xg = PhXq, we get 

Z\t) = Sh{T)X^+ f ShiT - s)PhF{X\s))ds 

Jo 

+ f Sh{T - s)Bh{X\s))dW{s), te [0,T]. (49) 

Jo 

Also let Z'^{t) = ShiT — t)X^(t), using fl44)) . we get 

Z'^it) = S,iT-t^)X:^ + {t-tJS,iT-t)Sl{t-t„,)P^FiX:i^) 

+ Sh{T - for t e [tmPm+l]. ^ ’ 


which can also be written as 


Z'^it) = Sh{T-tjX^+ / Sf,{T - s)Pf,F{Xi^)ds 

^ 'J tm 

+ f Sh{T-tjBh{X!^)dW{s), 

^ tm 


for i 6 


(51) 


Consequently, we obtain 

Z\T) = X^{t) = X^, ^"*(0) = .S;,(T)X'‘(0) = Sh{T)X^. 


Now the weak error is worked out as 

E[<h(X^)]-E[<h(X'*(T,Xo"))] = E[<l>(Z'*(T))]-E[/r^(T,Xo")] 

= E[u^{0, Z^{T))] - E[z/^(T, Z'*(0))] 

M-l 

= Y. - tm+1, Z’^{tm+l))] - E[z/'*(T - tm, Z^{tm))]- 

m=0 


Now using Ito formula (see Lemma [2.ip applied to G{t,x) = z/^(T — t, Z^{t)) in the interval 
and the fact that the ltd integral vanishes, we can write 


E[u>^{T - tm+l, Z^{tm+l))] - E[l^^{T - tm, Z>^itm))] 


^m+l 


E 


-iE 


tm+l 


{T-s,Z\s)) 


ds -\- 




E 


Dv\T - s, Z\s))Sh{T - s)PhF{Xt) 


Tr 


D^v\T - s, Z\s)) {ShiT - tm)BhiXt)] [ShiT - tm)BhiXY)Y 


ds 

ds. 


Using the fact that Du^it, ip) is identihed to an element of 14 (see the analogue representation 
at (HT]) 1 and (1451) . we hnally have 


12 









E[u^{T - tm+l, Z^{tm+l))] - E[z/^(T - tm, Z^{tm))] 


^m+1 


E 


'{Du^T - s, Z\s)), -5,(T - s)P,F(5,(s - T)Z%s)) + Sh{T - s)P,F(X^))h 


ds 


-|E 


^m +1 


Tr 


'dV(T - s, Z''(s))( {S4r - {SkiT - i„)Bk(X‘)}' 

-Si(T - s)B4A';i)(Sft(T - s)B^(At)y ] ds 


Pm+1 


E 


(P/(T - s, X'*(s)), P,P(X^) - P.P(X"(s))) 




ds 


-iE 


^m +1 


Tr 


P V(T - S, X\s)) ( {5,(s - trn)B,iX^)} {S,{s - t^)B,{Xt)y 


-B^{Xt)Bn{Xty 


ds. 


4. Weak convergence for a SPDE with additive noise 

In this section, we consider the additive noise where B = _ Jn order to prove our weak 

error estimate the following weak assumptions 0 o will be used. 

Assumption 4.1. [Assumption on nonlinear function F, and Q ] We assume that 


F : H ^ H is Lipschitz and twice continuously differentiable and satisfies 

||P(A)|| < L(1 + ||X||), XGP, (52) 

\\F'{Z){X)\\ < L\\Xl X,ZeH, (53) 

\\F"{z){x,,x,) < LIIX 1 IIIIX 2 II, Z, X, Xi, X 2 eP, (54) 

||(-A)-^P"(Z)(Xi,X 2 )|| < LIIX 1 IIIIX 2 II, Z,Xi,X 2 eP, /or some/3e [0,1), (55) 

||(-A)-ip'(Z)(X)|| < L(l + ||Z||i)||X||_i, (56) 

X e P, Z e Vii-Af/yfor some 7 G [1, 2). 

Furthermore, we assume that the covariance operator Q satisfies 

\\{-^)^Q^\\c 2 {h) < 00 , for some fi G (0,1]. (57) 


Note that the semigroup proprieties in Proposition 12.21 are satisfied for the discrete operator 
Ah. For our convergence proof, we add the following propriety to the discrete operator Ah. 


Proposition 4.1. Under Assumption |/. j| and Assumption 12. 11 the following proprieties are 
satisfied for discrete operators 


\\{-AhY-^^PhQ'^\\c 2 (H) < C+ Cph^ for some [3 (58) 

||(-A,)-^P,P"(Z)(Xi,X 2 )|| < C'(l + h2/')||Xi||||X2||, Z,Xi,X2GI4; /5g[0,1), (59) 
\\{-Ah)-"^PhF\Z){X)\\ < C'(P + (l + ||Z||i)||X||_i), Z,XeVh; ye [1,2), (60) 


^This assumption is weak compared to Assumption 12.2| and Assumption 12.31 
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Proof. According to (j^, Theorem 5.2) for all a > 0, the discrete operator and the 
continuous operator A satisfy 

\\Al^Ph-A-^\\LiH)<C^h^^, (61) 

where Ca is a constant dependent on a. For a = with (3 G (0,1], we have 

\\{A^-P,-A--)Q^\c,^h) = \\A^-PhQ-^-A-^Q^\c,^h) 

Now let be any orthonormal basis of H, since \\A~°'Q^\\c 2 (h) < oo, according to 

Assumption 14.11 frelation fl57)) h we have using fl6T]) 

\\Al-PnQ"^-A-<^Q^\c,iH) = E 

< WAh'^Ph- A~‘^\\li^h)\\Q"^\\c2{h) 


Now we can write, using flM]) and (ITT)) , that 

WA^PhQHciH) = ||AEnQ^-dl-“QUA-“Qi|U,(H) 
^ WAh'^PhQ^ — + ||A““Q 2 

< C + C^h^^ 

= C + Cah^-^. 


This proves (l58ll . Now to prove ([59]), we proceed as above and the fact that fl55ll is satisfied. 
Indeed for /3 G [0,1), using fIMl) . we have 

\\(Al^PH-A-f)F"{Z)(X,,X^)\\ < p^”P»-/l-'S||i,„,||||f"(Z)(A'i,A'2)|| < CLh^>\\XMX4 

Therefore, 

pyP,.F"(Z)(Ai, Ajjll < ||(AyP* - A-I>)F"(Z)(X,,X^)\\ + p-'’P"(Z)(Ai, A^jH 

< C(l + ft2'’)||.Y,||||.Y.||. 


The proof of fl60|l is done just like for fl5^ using fl56|l and fl5^ . ■ 

The estimates in Proposition 14.11 are very tight and can infiuence the order of convergence in 
space and time when jS in fl57|l is small. Indeed using the fact that 

h < C'(meas(f2))rf, (62) 


where meas(f2) is the either the length, the area or the volume of the domain 12 and d is the 
dimension, we have the following corollary. 


Corollary 4.1. Under Assumption\4- 1\ the following discrete proprieties are satisfied 


||(-A,)^p,g^lU,(H) < c, /3e(o,i]. 

||(-Ag-^P;,F"(Z)(Xi,X 2 )|| < CIIX 1 IIIIX 2 II, Z,Xi,X 2 GC;.; /?G [0,1), 
||(-A,)-ip,P'(Z)(X)|| < C(l + ||Z||i)||X||_i, Z,XGP,; ye [1,2), 


where the positive constant C is independent of h. 
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(63) 

(64) 

(65) 






The following lemma will be helpful for the proofs of convergence. 


Lemma 4.1. Assume that all the conditions above are fulfilled and that $ G C 2 (hf;M). For 
7 G [0,1], 7 i ,72 G [0,1) satisfying 7i + 72 < 1, there exists constants and c.^ 1,72 such that 

\\(-AhVDfi''{t,i>)\\ < (66) 

+ (67) 

where is defined by ([39j), and fj G I 4 . 

Proof. The proof is the same as in Note that although the linear operator in f^l is 

assumed to be self adjoint, the proof don’t make use of that. ■ 


The following results can be proven exactly along the line of 32|, Lemma 3.4 and Lemma 3.5]). 


Lemma 4.2. Suppose Assumption 12.21 is fulfilled and Xq G 'D((— then it holds for 
7 G [0, and arbitrary small e that 


sup \\{—AhyX^\\L.^(ji^H) < C and ||X^(f) — X'f^\\ 1 . 2 ( 0 ,h) < CAt 2 , 

0<m<M 

where X^f is defined by fimi and X^ is given in fl44D . Furthermore 

\\{-Ah)^Xlf\\L2(0,H) < <^(1 + At 2 ), 

for 13 G (0,1]. 

Proof. Before the proof of the lemma, we can prove exactly as in (sil . (3.3)] that 
ll^fclU 2 (B,H) < 00 , 0 < fc < M, sup \\PhF{Xj^)\\L 2 ( 0 ,H) < 00 . 

0<k<M 


( 68 ) 


(69) 


We concentrate on proving fl69]l . the proof of the other assertion can be done just as in 32 
Recall that 

+ A-p - I) P„F(XiW + (W',„ - . (70) 

(71) 


Bh = PhQ''"-Pk. 

This can also be written as 


rAt 


Xi = Si(A().Y;:._, + / S,(At - s)P,F{X;^Wds + S,(At)B,. . (72) 


'0 


Iterating gives 

Xt = ShitjX^ + / Snitn. - s)PHF{X'f)ds + ^,(f^ - h)BhAWk. 

1. n Ju 7. r. 


m —1 


(73) 


k=0 


Now 


-AA-2X’i 


M) < 


— Iq + h + h- 


m—1 r^k+l 

-A,,)-2Sh{trn)Xf,\\ + E / I 

k=0 J tp. 

m—1 

+ [^t Y1 \\{-Ah)^Sh{tm - tk)Bh\\l2lH) 

' k=0 


Ah)^Sh{tm — s)PhF{Xjj)\\L2{n,H)ds 
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( 74 ) 


We can prove exactly as in ( 32j, Lemma 3.5) that 

h < CAt^ and Jq < ||Wo||i < oo. 
Now it remain to prove that Ji is bounded above. We have 


m—1 

h<Y. \\{-Ah)--Sh{t^- s)PhF{Xl)ds\\L,i^^H) 

k={) Jtk 

^ C* Z] / {tm- s)~^{l + \\X]^\\L2{n,H))ds 
k=0 J t}^ 

< C (Yl f {tm-s)-^ds]( sup (1 + ||W^||L 2 (D,rr))) 

\k=0 Jtk / \0<A:<m—1 / 

= c([ sup (1 + 

\Jo / \0<fc<m—1 / 

< ( sup (1 + ||X^||i 2 (D,H))) < oo. 


Combining with fTMjl establishes the lemma. 


Now we can prove the convergence result for the semi disrete problem with non-self-adjoint 
operator A as the time step goes to zero. 

Theorem 4.1. Assume that Assumption 12.SI is fulfilled Xq G T>{{—Ay/‘^) and that $ G 
Then for arbitrary small e > 0, if the adjoint of the discrete operator A*^ satis¬ 
fies the analogue inequality as A^ in Corollary \f.l\ we have 

|E[<h(Wj^)] - E[<h(X^(T))] I < CAtf^-fi (75) 

where the constant C depends on fi, e, ...,T, L and the initial data, but is independent of h 
and M. 


Proof. We derived in fHHjl the following error representation 


M-l 


E[4.(.Y;,)| - E[4(A''‘(T))| = (6;„ + 6^) 


(76) 


m=0 


where the decompositions and 5^ are as follows 


i>L = 


tm+1 


E 


(Z1/(T - s, X^s)), P,F(X^) - FkF(Xys))} 


H 


ds 


(77) 


and 


5^ 


iE 


•'Tn+l 


Tr 


P V(T - s, X"(s)) ( - tm)Bh} - t^) - 


ds 


-|E 


^m+l 


Tr 


D^pi^T - s, Xys)){Sys - tm) - 


ds 


= 52,1 + 52,2 
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( 78 ) 













Note that our term 6^ is more simple than the one in Applying the Taylor’s formula to 
the drift term F we get 


F{X\s))-F{X^J = F\X^J{X\s)-X^J+ / F"{x{r)){X\s)-X^^,X\s)-X^J{l-r)dT, 

Jo 

where xi'^) = + r(X^(s) — X^), which allows to have 


\bln 


< 




E 


D^^\T - s,X\s)),P,F\Xt)iX\s) - X, 


H 


H. 


ds 


+ 




E 


\T-s,X\s)), / PnF"{x{r)){X\s)-Xt.X\s)-Xt){l-r)dr 


= F + F 

m ' m* 


H 


To estimate we use Corollary l4.1l Lemma 14.ll Lemma [4.21 and the Holder inequality 


m') 

^^m+l 


< 

^ m — 


E 




\S IX ..h 


-1 




E 


/ P,F'\xir))iX\s) - XI X\s) - X;^)(l - r)dr 
'0 ' ^ 
\^TX.,hfrTi vh/ 


ds 


-AhyDii\T-s,X\s)), 


{-Al)-^PnF"{x{r)){X\s) - X^ X\s) - X^)(l - r)dr 


H 


^^m+l 


< 


E 


{-AhfD^i\T-s,X\s)) 


X 


{-Al)-^P^F'\x{r)){X\s) - X:^, X\s) - X:y){l - r)dr 


h v"^ / 


ds 


ds 


< Cs 

< csL 


rtm+i 

I / E 

tm 0 

r^m+1 z*! 


E 


'tm Jo 

^^m+1 


{-Al)-^PnF"{x{r)){X\s) - Xy,X\s) - Xy) (T-syUrds 
X’^{s)-X!y ^]{T-s)-Urds 


< CXy-^ {T-s)-^ds 


We now turn to estimate Recall that from (jUj), 


(79) 


x\s)-xy = {sys-tj-i)xy + {s-tyy{s-tm)p,F{xy) 

+sys-tyByw{s)-w{ty). ^ ^ 


ds 
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Since the expectation of the Brownian motion vanishes, we therefore have 


71 < 




E 


(d/(T - s, X"(s)), PnF'iXt) {Sh{s - tm) - I)X^ 


H. 


ds 




E 


(T - s, X'^(5)), P^F\Xt){Sl{s - t^)PnF{Xt){s - P 


H. 


ds 


= I + IL 

Using Lemma [4.11 Corollary 14 .1 1 Lemma [4.21 and Proposition 14.11 yields 


(81) 




< C5 


E 


(-AD-lp^F'iX'J {S,{s - („.) - /).Y‘ (T - s)-i* 


rtm +1 r 


< c 


(i + lix;^lli)l||(5.(s-u-/)x;^ ^ {T-s)-lds 


r-Wl r 


< 


< 


1 + E||(-21,)-^)5,(5 - U - I){-A^)^Xt\\{T - sY^s 


r-tm+l r 


tn 


0 - 1 - 


1 + At— 


r^m+1 


1+0-e 


(s - E||(-A)—X^IKT - s)-i* 


< CAt"- / (T - s)-ids. 

Now we turn to approximate the term IF Note that we can rewrite 


as 


xYs)-xi = {sYs-tY-i)xt+ Shit - s)PhF{x';:^)dt 

^ tm 

+SYs -tYBYw{s) 


Therefore 
II = 

< 

< C 


pim+1 


E 


ptm+l rs 

I / E 

tm ^ tm 
ptm+1 /*-5 


XT-s,XYs)), / P.F'(X;^)5,(f - s)P,F(X, 


ds 


\\DYiT - s,X'‘(s))|| \\Pj,F\Xt)SYt - s)P,P(X, 


H} 

dtds 


E 


^ tm ^ tm 

< CE[\\PhF{X, 

< CA0. 


\\SYt-s)P^,F{X, 

rtm+l ps 

I] / / dtds 

tm tm 


dtds 


Combining fITB]) . fl52D and fl52]) . we have 

|6;„| < CAi«- 


P^m+l 


(T - sY^ds + CAf 

M-l 

\bln\ < CAY-^ 

m=0 
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(82) 


(83) 


(84) 








































The term 6^ can be approximated just as in [32| to be 


M-l 

Y, \bl\ < CAtO-. 

m=0 


(85) 


Finally we therefore have 

|E[<F(Xj^)] - E[<F(X^(T))] I < (86) 


5. Weak convergence for a SPDE with multiplicative noise 

Now we consider the case of a stochastic partial differential equation with multiplicative noise, 
that is ([I]) 

For convergence proof, we make the following assumption on the noise term. 

Assumption 5.1. We assume that there exists a constant a > 0 such that 

\\{-A)^B{X)\\c,(H)<C{l + \\X\\p.,), forallXeH, /3e[0,l]. (87) 

One can prove just as in the proof of Proposition 14.11 that the discrete operators A^, and 
PhB{X) satisfies 


\\{-Aj,)^PhB{X)\\c,(H) < Ch^-^il + ||X|l;3-i) < C"(l + ||X|l;3-i), for all X e H. (88) 
For (3 G [0,1], we also have using Proposition 12. II 


\\{-Ah)~ Bh{X)\\c2(H) = 

< 

< 

< 


||(-A)^P,5(X)P,|U,(^,) 
C\\{-A^)^P^B{X)\\c,(h) 
O'(l + ||X||^_0, 

C"(1 + ||X||), forallXGP. 


(89) 


Thanks to Lemma |T]2] and fl88l) (or fl8^ 1 all the results presented above for additive noise also 
hold for multiplicative noise and we then have the following convergence result. 

Theorem 5.1. Assume that all the conditions of Theorem \4.1\ are satisfied with the condition 
fl35]l replaced by fISS]) . For Xq G P((—A)^/^), $ G C 2 (P;M) and for arbitrary small e > 0, i/ 

the adjoint of the discrete operator Al^ satisfies the analogue inequality as Ah in Corollary \j.l\ 
we have the following weak error convergence rate 


|E|4(A:!i)] - E|4(X'“(T))]| < CAC-, 


(90) 


where the constant C depends on fi, e, ...,T, L and the initial data, but is independent of h 
and M. 
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Proof. The proof follows the same lines as that of Theorem 14.11 


M-l 


E|4(A';,)| - E|4(A'*(T))| = (t + 


( 91 ) 


m=0 

The term 6^ is the same as in Theorem 14.11 Here 6^ is given by 


6^ = lE 
2 




Tr 


v/ tn 


+ |E 

^ tm 

= 52,1 ^ 52,2 

m ' m 


^^m+l 


- s, X^{s)) {Shis - tm)BhiX!^)) 

x{iShis-trr,)-I)BhiX^)y]ds 

- s,Xys))iShis -tm) - I) (BhiXt)) (5/i(X^))^ 


ds 


(92) 


Let us estimate b‘^. Proposition 12.11 allows to have 
b^j} 

D V(T - s, [Shis - tm)BhiX^)] 

X [iShis-t^)-I)BhiX:i^)]'- 


< iE 


v/ tn 


= iE 


P^m+1 


i-Ah)^-^D^fi\T - s,X\s))i-Ah)^i-Ah)^ [Shis-tjBhiX^)] 


CiiH) 

/ 3-1 


ds 


X 


i-Ah)-^^%Shis - tm) - I)BhiX^) 

Using Lemma im Proposition l2Tl Proposition lT^ and fl8^ yields the following estimations 


- !(C 

ds 

- 

Ci(H) J 


< (PE 


< CE 

X 

< (PE 


^^m+1 




13+1 


i-Ah)^Shis - t^)BhiXt)\i-Ah)^-%Shis - tj - I)BhiX, 

x[iT - s)-^+^ + l]ds 


Ci{H) 


-Ah)^Shis-t^)BhiXl^) 


Cc(H) 


(-AJ-'¥+.(S„(s - („) - /)B»(.Y‘) - s)-'+‘ + l]<is 

tm+1 g_^ 

-^'•)^^'-(^™)IL„„, 

X (-Ai)-'’+‘(Sk(s - «„) - |(T - s)-‘+' + l]£is 

C2iH) 

P^m+1 

< C / sup E(l + ||Xf 11)2. [(T-s)-'+^ + l]ds 

Jtm 0<i<M L{H) 

r^m+i 

< iT - s)-^^^ds + CA+-^+\ 


20 
































































Similar as for we have 

^2,2 

rrt 


= 


trn+l 


< CE 


Tr 

-i+l 


x(-.lj)-“f"+-(S^(s - t„) - /) ((-A)'¥Bk(A‘) 


(is 


x[(T-s)-^+" + l](is 

< / "^\T-s)-^+^rfs + CAt^-^+^ 


^m+1 


£i(H) 


So by summing up as in Theorem 14.11 the proof is ended. 


6. Strong convergence and toward full weak convergence results 

The goal here is to provide the space and time convergence proof of the exponential scheme. 
Before that we will provide strong convergence of the semi discrete solution. 

Theorem 6.1. Let and Xh be the solutions respectively of fl35l) and the following semi 
discrete problem 


dXh = {AhXh + PhF{Xh))dt + PhB{Xh)dW (93) 

A;,(0) = PhXo. 

Let {3 G [0,2), assume that Assumption \2fM and Assumption \2.S[ are satisfied. For (3 G 
[0,1] assume that the relation fl57|l of Assumption \4^ (when dealing with additive noise) 
and Assumption 15. il (when dealing with multiplicative noise) are satisfied. For (3 G [1,2) 

assume that B{T>{{—A)^)) C HS (^^F{H),T>{{—A)^'^ and ||(—y4)'^i?(n)||j ;^2 < c(l + 

||n||^_i) for V G V{{—A)^). If Xq G L 2 {E),'D{{—A)^F)), there exists a positive constant C 
independent of h such that the following estimations hold: 

\\X{t)-XhmL,(p,H) < (3 e [0,1]. (94) 

Furthermore assume that the linear operator A is self adjoint, the following estimation hold 

\\Xit) - X\t)\\L,(n,H) < Ch^, /S e [0,2). (95) 


Before prove Theorem 16.11 let us make some remarks and provide some preparatory results. 


Remark 6.1. Theorem 16'. il extends |7i. Theorem 1.1] for non-self-adjoint operator A and 
also provide optimal convergence proof for (3 G [0,1) which was not studied in li.W . 
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Remark 6.2. For additive noise, we can observe from 113] that the semi discrete problem 
fl35|) is equivalent to the following problem, find Xh{t) G Vh such that 

dXh = {AhX^+ PhF{X^))dt + dWh (96) 

X'^(O) = PhXo. 

where Whit) is a PhQPh-Wiener process on Vh with the following representation 

Nh 

Wh{t) = y^^y/^eh,i(di{t), (97) 

i=l 

where {qh,i,eh,i) are the eigenvalues and eigenfunctions of the covariance operator Qh : = 
PhQPh and Pi are independent and identically distributed standard Brownian motions. More 
precisely {qh,i,eh,i) is the finite element solution of the eigenvalue problem Qu = 7 W. 

If the exact eigenvalues and eigenfunctions (g*, ep of the covariance operator Q are known, 
replacing in fl96]) (or fl5Sl) ) Wh{t) by Pk^'‘(t), defined by 

Nh 

'VH) = JIViteMt), (98) 

i=l 

will not necessarily change the optimal convergence order in our scheme. From fT^l . it is 
also proved that if the kernel of the covariance function Q is regular and the mesh family is 
quasi-uniform, it is enough to take M < Nh noise terms in fl98D (or ra; without loss the 
optimal convergence order. 

Of course for multiplicative noise PhW can also be expanded on the basis ofVh with Nh terms 
(see fM]). 

As we are also dealing in Theorem 16.11 with non-self-adjoint operator in fl9T|l . let ns provide 
some preparatory resnlts before giving the proof of Theorem 16.11 

We introdnce the Riesz representation operator Rh '. V ^ Vh dehned by 

{-ARhV, x) = {-Av, x) = a{v, x), n G R, Vy G I 4 . 

Under the usual regularity assumptions on the triangulation and in view of U—ellipticity fl2T]l . 
it is well known (see [sj) that the following error bounds holds 

\\RhV - v\\ + h\\RhV - v\\Hi(n < CN\\v\\Hr-{n), veVnH^{Il), rG{l,2}. (99) 

By interpolation, we have 

WRhV — v\\h\\RhV — < CN\\v\\r, vE'D(A''P), 1 < r < 2. (100) 

Let us consider the following deterministic problem, which consists of hnding u E V such that 
such that 


u' = Au given m(0) = v, 
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fG (0,T]. 


( 101 ) 



The corresponding semi-discretization in space is : Find Uh G i4 such that 

Ufi AfiUfi 

where = PhV. Dehne the operator 

nit) := Sit) - Shit)Ph = - e^^'^Ph ( 102 ) 

so that M(t) — Uhit) = Thit)v. 

The following lemma will be important in our proof. 

Lemma 6.1. The following estimates hold on the semi-discrete approximation of fllOip . There 
exists a constant C > 0 such that 


(%) Forve V{{-A)^/^) 

\\u{t) — Uh{t)\\ = ||T/i(f)n|| < l<r<2, 0 < 7 <r. (103) 

(%t) Forve T>((-A)('^-0/2) 


(*t \ 2 

|2 


\\This)v\fds ) < C'h'’'||t>||..y_i, 0 < 7 < 2. 


(104) 


Proof. The proof of (i) can be found in j22| using fllOOp . For self adjoint operator, the proof 
of (ii) is done as in Lemme 4.1] if 7 G [0,1] and in[l3, Lemma 4.2] if 7 G [1,2) and the 

parameter r used in [l9|. Lemma 4.2] is r = 7 — 1. Both proofs only uses general concepts and 
not spectral decomposition of the linear operator A, so can easily be generalized. Note that 
the non-self adjoint case should make use of M, 4.17] or j^. Lemma 4.3] instead of 31, 2.29] 
used in [l^. Lemma 4.2]. 


Let us now provide the proof of Theorem 16.1 


Proof. The proof of the estimation fl95|) for additive noise can be found in 1^ and can be 
updated to multiplicative noise following (jjj ]. 

For /3 G [1,2), the proof for the estimation fl9T|l when A is self adjoint operator can be found in 
flij ] as the Assumption 12.21 implies jl^. Assumption 2.1] by taking r = — 1 in [l^. Theorem 

1.1]. Let us give more general proof by closely follow 18|, ll9|. The corresponding mild solution 
of (1951) is given by 


Xnit) = Shit)PhXo+ / Shit-s)PhFiXhis))ds+ / Shit-s)PhBiX^,is))dWis). (105) 

Jo Jo 

Indeed, we have 

||X(f) - Xhit)\\L2(0,H) 

< \\mt)PhXo+ f Sit-s)FiX{s))ds- f Shit-s)PHF(Xhis))ds\\L 2 iB,H) 


+ 11 / Sit - s)BiXis))dWis) 

Jo 

= h + h. 


Shit - s)PhBiXhis))dWis)\\L,(^^,H) 
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(106) 















The estimation of Ji is the same as in 1^ and we have 


Ii<C ||X(s) - Xh{s)\\L,iB,H)ds + Ch^. 
Jo 


For the estimation of I 2 , we follow closely [19j. Indeed we have 


\\S{t - s)BiX{s)) - S,it - s)PhB{X,,{s))\\l^^^)ds 


I2 < (FIE 
< C 

+c 

+c 


1/2 


\\S,it - s)P,iBiX,is)) - BiXismi^^j^^ds] 


\ 1/2 


'0 


\\n{t-s){B{X{s))-B{Xm\\l^^)ds 


1/2 


Z/2 (ID),M) 


'■* \ 1/2 

11^(4-s)B(A'(«)|||,,„,* 


Z/2 (OjM) 


= n + n + ii 


(107) 

(108) 

(109) 

( 110 ) 

( 111 ) 

( 112 ) 


The stability propriety of the semi group Proposition 12.21 and the Lipschitz condition in As¬ 
sumption 12.21 allow to have 


(113) 


l2<C[ l|A(s)-A4s)||i,,„_„|* 


Following closely [19|, but with fllOdp in Lemma [ 6 . II with r = /3, 7 = 0, allow to have 

/| < Ch^. (114) 


For /3 G [0,2), as in 19|, by using (ii) in Lemma [6.11 gives 

< ChJ. 

Coming back to / 2 , we have 

rt rt 


(115) 


10 

,2-e 


< Ch'^-^ + C{ I \\X{s)-X,{s)\\l^^^^)ds 


< Chf^ + C{ I \\X{s)-X,{s)\\l^^^^^ds 


S{t - s)B{X{s))dW{s) - / Sh{t-s)PHB{Xh{s))dW{s)\\L,io,H) (116) 

Jo 

'' ’ (117) 

(118) 

Combining Ji and I 2 gives 

||A(i)-A^(«)||7„_„, < Ctf^ + C f\\X(s)-Xu(s)\\l,^„_„/s. (119) 

Jo 

Gronwall’s lemma is therefore applied to end the proof. ■ 

The following theorem provide the full weak convergence when the solution is regular enough. 
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Theorem 6.2. Let X and be respectively the solution of ([T]) and the numerical solution 
from fl38|) at the final time T. Let fd G [1, 2), assume that Assumption \2.^ (for multiplicative 
noise, all conditions except m) and Assumption IS.A (for multiplicative noise) are satisfied. 
For fi = I (trace class noise) assume fl57|) of Assumption M (when dealing with additive 
noise) and Assumption \5 . 1\ (when dealing with multiplicative noise) are satisfied. For fd G (1, 2) 
assume that fl57|) of Assumption \4A\ is also satisfied for additive noise, but B{V{{—A)^~)) C 
HS and, \\{-A)^B{v)\\l 2 < c(l + HvH/j-i) for v G V{{-A)^) 

for multiplicative noise. Furthermore assume that Xq G P ((—and the linear operator 
is selfadjoint. For $ G and arbitrary small e > 0, the following estimation hold 

|E[<I>(X]^)] - E[$(X(T)]| < + h^), (120) 

where the constant C depends on a, e, L and the initial data, but is independent of h 

and M. 


Proof. Indeed we have the following decomposition 


\E[^X(Ij)] -E[<I)(X(T)]| < |E[$(X)(^)] -E[<I>(X'^(T)]| + |E[<I>(X'*(T)] - E[d>(X(T)]|. (121) 

Note that if Xq G 'D{(—A)h/^), 1 < (d < 2, then Xq G 'D((—A)^/^) and all the lemma nsed in 
Theorem 14.11 and Theorem 15.II are still valid. From Theorem 14 . 1 1 and Theorem 15.11 with /d = 1 
(as the noise is assnmed to be trace class), we have 


\E[^{X’f^)] - E[<I)(X^(T)]| < CAt 


l-e 


( 122 ) 
3. Re- 


Note that this temporal order is the double of the strong convergence obtained in 
member that we are using low order hnite element method for space discretization, where the 
optimal convergence is 2 (for deterministic case), therefore for Xq G T>{{—A)^(^), 1 < (d < 2, 
according to Theorem 16.11 and Remark 16.21 we cannot expect the order of weak convergence 
to double the strong order in space. The weak convergence order in space (of course, not 
necessarily optimal) can be the same as the strong convergence order when the solution is 
regular enough. Using the fact that <I> G so is Lipschitz, we therefore have 


|E[<I)(X'*(T)] -E[<I)(X(T)]| 


< |E[$(X'*(T)) - <I>(X(T))]| 

< E||X^(T)-X(T)|| 

< C'||X'^(T)-X(T)|U,(„,h) 

< Ch^. 


The following theorem using recent result in the literature shows that the space order of 
convergence in Theorem 16.21 is far to be optimal for fd = 1. 

Theorem 6.3. Assume that A = A with Dirichlet boundary condition (V = Fl^lfiV)), assume 
that the noise is additive, Assumption \f.l\ is satisfied with \\{—A)^Q^^‘^\\c 2 {h) < oo for some 
(d G [1/2,1]. Furthermore assume that Xq G V{{—A)^). For $ G C 2 {H;E) and arbitrary small 
e > 0, the following estimation hold 

|E[<h(X](f)] - E[<h(X(T)]| < C{At^-^ + 
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(123) 















where the constant C depends on a, e, L and the initial data, but is independent of h 

and M. 

Let us prove Theorem 16.31 

Proof. The proof follows the one for Theorem 16.21 and we have 

|E[<h(XjJ:^)] - E[<h(X(T)]| < |E[<h(Xj^)] - E[$(X'‘(T)]| + |E[<h(X^(T)] - E[<h(X(T)]|. (124) 

If Xq G V{{—A)h), 1/2 < (3 < I, then Xq G P((—and all the lemma used in Theo¬ 
rem o and Theorem 15.11 are still valid. From Theorem 14.11 if ||(— Q^/‘^\\c 2 (h) < oo for 
some (3 G [1/2,1] we have 

|E[<F(Xj(:^)] - E[<1)(X^(T)]| < CXt^-\ (125) 

Note that Assumption 14.11 implies that F G C\{H]H). For \\{—A)^Q^/‘^\\c 2 {h) < oo for some 
/3 g [1/2,1], B Assumption A, Theorem 1.1] gives 

|E[$(X^(T)] - E[<F(X(T)]| < (126) 

The proof of (11261) in uses some elements of Malliavin calculus. ■ 

The following remark generalizes the Theorem 16.31 for general selfadjoint operator with not 
necessarily Dirichet boundary condition. 

Remark 6.3. For additive noise and under the same condition as in Theorem, \6.A if A is 
selfadjoint, and \\{—A)^Q^/‘^\\c 2 {h) < oo for some f3 G [k, 1], k > 1/2 and Xq G V{{—A)h). 
For $ G and arbitrary small e > 0, the following estimation hold 

|E[<F(Xj(^)] - E[<F(X(T)] I < C{At^-^ + . (127) 

The proof is the same as in Assumption A, Theorem 1.1] where the Laplace operator is 
used just for simplicity. The difference comes from the set of the eigenvalues of the self adjoint 
operator A. Once the range of (3 such that || (—A)^Q^/^||£ 2 (rr) < oo is found, the proof ofU> 
Assumption A, Theorem 1.1] is applied line by line. 

The following numerical simulations will conhrm numerically estimation fll27p of Remark 16.31 


7. Numerical Simulations 


We consider the reaction diffusion equation 


dX = {DAX - 0.5X)dt + dW given X(0) = Xq, (128) 


on the time interval [0, T] and homogeneous Neumann boundary conditions on the domain = 
[0, Li] X [0, L 2 ]. The noise is represented by ([2]). The eigenfunctions {eij}ij>o = 
of the operator —A here are given by 





cos(Af^x), 


A® = 0, 



ill 

u 
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where / G {1, 2}, a; G and i G N with the corresponding eigenvalues given by 

Aij = (aP)^ + We take Li = L 2 = 1. Notice that A = DA does not satisfy 

Assumption 12.11 as 0 is an eigenvalue. To eliminate the eigenvalue 0 we use the perturbed 
operator A = DA + el, e > 0. The exact solution of fll28|) is known . Indeed the decomposition 
of 01281) in each eigenvector node yields the following Ornstein-Uhlenbeck process 

dXi = —{DXi + 0.5)Xj(if + y/^idl3i{t) i G (129) 

This is a Gaussian process with the mild solution 

W(t) = e-'='%(0) + ki = DXi + 0.5, (130) 

Jo 

which is therefore an Ornstein-Uhlenbeck process. Applying the Ito isometry yields the fol¬ 
lowing exact variance of Xi{t) 

Var(W«)) = ^(l-e-“-‘). (131) 

During the simulations, we compute the exact solution recurrently as 

pt'm + l 

xp+l ^ ^-k,Atxjn + ^ / 

/ \ 1/2 

where Ri^m are independent, standard normally distributed random variables with mean 0 
and variance 1. The expression in 0132p allows to use the same set of random numbers for 
both the exact and the numerical solutions. 

Our function F{u) = —0.5m is linear and obviously satishes the Lipschitz condition in As¬ 
sumption 12.21 We assume that the covariance operator Q and A have the same eigenfunctions 
and we take the eigenvalues of the covariance operator to be 

= (133) 

in the representation ([2]) for some small h > 0. Indeed for (3 G [0,1] we have 

||(-A)VqV 2||^^{^) < oo ^ < TT^ ^ if + < co. 

hJ)6N2 (iJ)GN2 

We will consider the following two test functions 

: / ^ / f{x)dx, $2 : / ^ ll/ir. (134) 

Jn 

which obviously belong to By setting B = DA — 0.5/, using the fact the the Ito’s 
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integral vanishes, for Xq G if we have 


E$i(X(i)) = / e^^Xo{x)dx + E 


e^t-^^^dW{s,x)) dx 


n \Jo 


Xo(x)dx + / E 


'o L~'o 


„(t-s 


^^dW{s, 


X) 


dx 


Xo{x)dx 


in 

OO 


/ eij{x)dx, kij = DXij + 0.5. 

i,j>0 


We also have 

E42(A'(«)) = E||e‘®.Yo+ f 


( 135 ) 


(136) 


= E||e‘®.Y„||'^ + E|| / e(‘-*)®diy(i>)f+ 2E{e‘®.Yo, / e(‘-*)®diy(i>)),j(137) 
Jo Jo 

Using the spectral decomposition, the fact that Xq E H and the Ito isometry yields 


OO 

Ed)2(X(f)) = — (e 

i,j>0 


i,j>0 


H- 


(138) 


To compute the exact solutions E<hi(X(i)) and Ed> 2 (X(t)), we can either truncate fll35p and 
fll38p by using the hrst Nh terms, Nh being the number of hnite element test basis functions, 
or use directly the Monte Carlo method with fll32p . 

Our code was implemented in Matlab 8.1. We use two different intial solutions with each test 
function <hi and <h 2 - Indeed we use the initial solution Xq = Xq = 0 for $ 2 , and 


.Y„ = .Y<‘' = 

i,j>l 




-1.001 


for for $ 1 . Our hnite element triangulation has been done from rectangular grid of maximal 
length h. Fast Leja Points technique as presented in j^, 1^ is used to compute the exponential 
matrix function = (fi. In the legends of our graphs, we use the following notations 

• ’’TimeErrorl” denotes the weak error for hxed h = 1/150 with <F = $ 1 , Xq = Xq^^ and 
hnal time T = 1. 

• ”TimeError2” denotes the weak error for hxed h = 1/50 with $ = $ 2 , = -^ 0 ^^ and 

hnal time T = 1. 

• ’’SpaceErrorl” denotes the weak error for hxed time step At = 1/500 with $ = $ 1 , 
Xq = Xg^^ and hnal time T = 0.1. 

• ” SpaceError2” denotes the weak error for hxed time step At = 1/20000 with $ = $2 
and Xg = Xg^^ and hnal time T = 0.1. 












Figure 1: (a) Weak convergence in time of the exponential scheme for $ = $i with Xq = (TimeErrorl) 
and <& = <1)2 with Xq = Xq ^ (TimeError2). (b) Weak convergence in space of the exponential scheme for 
<1) = (hi with Xq = Xg^^ (SpaceErrorl), and <1) = <1)2 with Xg = Xg^^ (SpaceError2). 


In all our graphs D = 0.1, /3 = 1, 6 = 0.001, the weak errors are computed at the hnal time T 
and 50 realizations are used to estimate the weak errors with the Monte Carlo method. 


Figure 1(a) shows the weak convergence in time of the exponential scheme for <I) = $i with 
Xq = and <h = <1)2 with Xq = Xq^\ The weak order of convergence in time is respectively 
1.1 for <h = $1 and 1.008 for <I) = <1)2. These orders of convergence are then closed to the 
optimal order 1 obtained in Remark 16.31 


For space convergence, very small time steps are needed and the weak errors are performed at 
small hnal time T = 0.1, Figure 1(b) shows the weak convergence in space of the exponential 
scheme for <!> = <I)i with Xq = Xq^\ and <!> = <1)2 with Xq = Xq^\ The weak order of 
convergence in space is respectively 1.71 for $ = <I)i and 2.01 for <h = <h 2 . These orders of 
convergence are also closed to the optimal order 2 obtained in Remark 16.31 To sum up, our 
simulations conhrm the theoretical results obtained in Remark 16.31 
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